Evaluation of Shear Strength and Stiffness of a Loess–Sand Mixture in Triaxial and Unconfined Compression Tests

Mechanical soil parameters are not constants and can be defined in various ways. Therefore, determination of their values for engineering practice is difficult. This problem is discussed based on results of piezoceramic element tests and triaxial tests (unconfined and confined) on loess specimens improved by compaction and sand admixture (20% by weight). The study indicated also the effectiveness of this simple method of loess stabilization. The influence of specimen size, draining conditions, stress and strain state, and different calculation methods on the evaluation of basic mechanical parameters were analyzed. The initial shear and Young’s moduli, the degradation of secant moduli with strain, tangent moduli, and Poisson’ ratio were determined. The results showed that the shear strength parameters are much less sensitive to the test variables than the stiffness parameters are. In triaxial tests, the strength criterion adopted, the sample size, and the drainage conditions influenced the measured value of cohesion, with a much smaller impact on the angle of internal friction. On the other hand, the adopted definition of the parameter and the range of strains had the greatest influence on the value of the stiffness modulus. Moreover, larger specimens were usually found to be stiffer.


Introduction
Evaluation of soil shear strength and stiffness is one of the main tasks of geotechnical investigation.These properties are crucial at every stage of the 'lives' of geoengineering structures for their safe and economical design, execution, monitoring, improvement, or demolition.Due to the complicated nature of soils, the mechanical parameters cannot be treated as constants [1], as they depend on several factors, such as current stress and strain conditions, the genesis of the deposit, mineralogical composition, loading history and loading velocity, saturation degree, etc.Even if the study is limited to the results of laboratory tests, representing just one point of the subsoil, the data can be interpreted in various ways, and, thus, give different end results.As there are no well-defined common guidelines for the selection of the particular test methods and calculation procedures, this may be very confusing for civil engineers.
The required set of parameters needed to describe the complete soil behavior depends on the constitutive model that is to be used in the design [2].However, the basic set, found in practically every technical document and used in most models, includes the following parameters: cohesion c and angle of internal friction φ to define the shear strength boundary, and shear modulus G or Young's modulus E, and Poisson's ratio ν, for soil stiffness in the pre-failure stage.They all can be evaluated in basic field or laboratory tests, although the definitions of these parameters and the conditions in which they should be determined are not always clear.
In this paper, the parameters of shear strength and stiffness, i.e., c', φ', E, G, and ν, were evaluated, based on results of conventional triaxial and unconfined compression experiments, supplemented by piezoceramic element tests.The study was carried out on a saturated, remolded, and compacted loess with 20% sand admixture (by weight).This particular soil was chosen because the pre-failure behavior and shear strength of this material have not been well recognized yet.Compared to the research works on undisturbed loess [3][4][5][6][7][8][9], there are much fewer studies referring to the mechanical parameters of the remolded and compacted loess [7,[10][11][12][13][14][15].Several studies focused on chemical stabilization, e.g., [16][17][18][19], but very few papers [20,21] refer to loess improved by the addition of a coarser fraction (e.g., sand), which is one of the simplest, yet poorly recognized, solutions.Researchers usually concentrate on selected types of behavior and parameters.For instance, (a) Zhang et al. [19], Ma and Ma [22], and Li et al. [23] investigated the unconfined compression strength of compacted loess; (b) Kim and Kang [11] and Zhang et al. [19] obtained direct shear strength parameters; (c) Zhong and Liu [9], as well as Li et al., [23] tested undisturbed loess in triaxial apparatus to determine its behavior and shear strength; (d) Ma et al. [12] and Xu et al. [15] assessed c' and φ' in triaxial tests on consolidated specimens that were sheared in undrained conditions; (e) Rinaldi et al. [5,6] and Wang et al. [8,13] investigated the small-strain shear modulus G 0 under isotropic consolidation, using bender elements; (f) Liu et al. [20] determined the effect of water and sand content on small-strain stiffness, using a resonant column; (g) Song et al. [7] examined the G 0 and damping of undisturbed loess in torsional resonant column tests; (h) Bowers [24] evaluated Poisson's ratio and the initial tangent Young's modulus for unsaturated loess.The literature results point out the influence of grading, water content, and void ratio on the material behavior.Contrary to the mentioned analyses, in this study, special attention was paid to the various factors influencing the soil strength and stiffness, such as saturation (low or complete), confining stress (50-400 kPa), strain (ca.0.01-10%), draining conditions, specimen size, and strain measurement method.This allowed a comprehensive comparison of the mechanical parameters evaluated based on different test procedures and interpretation techniques.Such an approach is rare, yet it is necessary for better understanding of material behavior [25].In addition, the majority of studies refer to the loess from the Chinese Loess Plateau [7][8][9][12][13][14][15]19,20,22,23]; only few consider loess from America [5,6,11,24], and much less is known about the European loess.Therefore, the results obtained in this research for loess deposited in Poland are filling a knowledge gap.Some preliminary outcomes referring to the stiffness of the loess-sand mixture have been presented in [21].This study is an extension and continuation of that work.The results are also compared with those published for similar soils.

Mechanical Parameters Interpretation
Cohesion and the angle of internal friction are the parameters (initial value and slope) of a linear Coulomb function, approximating the failure boundary in the shear stressnormal stress (τ-σ') space.The boundary is generally curvilinear [26], thus the linear approximation is valid only for a given range of σ.The failure boundary can be defined in various ways-it can describe, e.g., the peak shear strength, as well as critical (steady), unstable, or residual states, giving different values of the parameters.This also means that the cohesion and angle of internal friction should not be treated as intrinsic soil properties.It is worth mentioning that c = 0 kPa, according to the critical state theory, which is also consistent with the approach presented by Santamarina [27] and Schofield [28].
Although commonly used in all design situations, Young's modulus and Poisson's ratio are, theoretically, relevant only there, where the minor σ 3 and intermediate σ 2 principal stress tensor components are zero.This can be simulated in unconfined compression tests.The parameters are calculated with the following Equations ( 1) and ( 2): where ∆σ' 1 and ∆ε 1 are, respectively, the increments of axial (vertical, major) effective stress and strain, and ∆ε 3 is the increment of strain in the direction perpendicular to the loading (horizontal, minor).Unconfined compression rarely occurs in a real subsoil subjected to loading and is simply infeasible for some soils (e.g., loose sands).Therefore, in geotechnics, triaxial tests on confined specimens are much more preferred.In triaxial compression, the horizontal components of the strain and stress tensor are equal, so, based on the theory of elasticity [29], the Poisson's ratio and the Young's modulus during shearing, at ∆σ' 1 > ∆σ' 3 , shall be estimated with the following Equations ( 3) and ( 4): where ∆ε 1 , ∆ε 2 , ∆ε 3 and ∆σ' 1 , ∆σ' 2 , ∆σ' 3 are the major, intermediate, and minor principal strains and stresses, respectively.During an undrained shearing, the Equation ( 4) reduces to Equation (1), because the Poisson's ratio ν can be assumed to be equal to 0.5, while in a drained test, as ∆σ' 3 = 0, Equation (4) reduces to the following Equation ( 5): where ∆q is the increment of effective or total stress intensity, being the difference between the major and minor principal stress (q = q' = σ' 1 − σ' 3 ).The E values depend on the drainage conditions, so two types of elastic moduli are usually distinguished-'drained' E' and 'undrained' E u .
The shear modulus G is theoretically defined as the ratio of shear stress τ and shear strain γ in the conditions of simple shearing, i.e., with zero normal stress tensor components (mean effective stress p' = 0) and with no volume change (ε vol = 0).It is considered independent of drainage conditions [30].Although the true simple shear can be simulated only in a simple shear apparatus [31], G is often evaluated in triaxial compression tests, using the following Equation ( 6): where ε s is the strain intensity calculated as 2/3 of the difference between the major (vertical) ε 1 and minor (horizontal) ε 3 principal strains.The influence of volumetric strain can be treated as negligible in undrained conditions.In such a case, ε s corresponds to ε 1 and G to 1/3 of the elastic 'undrained' modulus E u .However, in a standard test with constant confining effective stress, this approach does not meet the remaining simple shear condition of ∆p' = 0.Both the shear and elastic moduli can be calculated as the secant or tangent (c.f. Figure 1 in [21]).They strongly depend on the current values of stress and strain, and on the previous loading history.In triaxial compression, the secant modulus is always positive, while the tangent modulus decreases to zero at the maximum stress intensity q f and becomes negative if the soil softens with strain.Constitutive modelling often uses a specific value of the secant modulus obtained in a monotonic loading at 50% of q f , named E 50 .
stratigraphy can be found in [38][39][40].In this area, the loess deposits are typically 4-6 m thick and contain large amounts of sand and clay fractions, with organic matter and calcium carbonate contents of approximately 1% and 0-5%, respectively.The soil sample was taken directly from a slope exposure.From the microscopic image shown in Figure 1, it is clear that the silty particles are typical for loess-they are mostly isometric and consist of quartz crystal fragments.The material was passed through a 1 mm sieve, dried, and mixed with a medium quartz sand in the proportion of 80%:20% by weight.The grain size distribution curves of It is well known that the soil stiffness decreases with strain [32,33] and that various strain measurement techniques have to be used to accurately describe the degradation in triaxial tests [21,25,34].A complete set includes external and local (on-sample) strain transducers to read large (>10 −3 ) and small (10 −6 -10 −2 ) strains, respectively, and a pair of piezoceramic elements (bender-extenders) for the very small strain range (<10 −5 ).Piezoceramic elements (PE) allow non-destructive measurement of the velocity of shear wave (S-wave) and compression wave (P-wave), V S and V P , passing through a soil specimen.If the current bulk density of the soil ρ is known, the initial (maximum) shear modulus G 0 can be calculated based on the following Equation ( 7): It is believed that PE generate strains smaller than 0.001% [35] and that in this range soil behaves as a linearly elastic material.Based on the theory of elasticity and assuming that the soil is isotropic, the initial Young's modulus E 0 can be calculated from the following Equation ( 8): There is no distinction between the secant and tangent initial moduli.The initial Poisson's ratio ν 0 can be determined from the PE tests [36] based on the following Equation ( 9): It shall be noted that the P-wave velocity depends strongly on the degree of saturationtherefore, this method of evaluation of the Poisson's ratio is not recommended for saturated soils.

Materials and Sample Preparation
The study was carried out on a loess from the Sudety forelands deposit in Trzebnica Hills, in the neighborhood of the city of Wrocław (Poland).The site belongs to the northern European loess belt [37].Studies on its sedimentology, paleopedology, and stratigraphy can be found in [38][39][40].In this area, the loess deposits are typically 4-6 m thick and contain large amounts of sand and clay fractions, with organic matter and calcium carbonate contents of approximately 1% and 0-5%, respectively.The soil sample was taken directly from a slope exposure.From the microscopic image shown in Figure 1, it is clear that the silty particles are typical for loess-they are mostly isometric and consist of quartz crystal fragments.
The material was passed through a 1 mm sieve, dried, and mixed with a medium quartz sand in the proportion of 80%:20% by weight.The grain size distribution curves of the sand, loess, and the sand-loess mixture (after the removal of calcium carbonate; ISO 17892-4 [41]) are shown in Figure 2. The physical parameters of the mixture, namely sand (Sa), silt (Si), and clay (Cl) contents, mean grain size d 50 , coefficients of uniformity C U and curvature C C , plastic limit w P , liquid limit w L [42], plasticity index I P = w L − w P , colloidal activity A, loss on ignition LOI at 600 • C [43], particle density ρ s [44], optimum moisture content w opt , and maximum dry density ρ d .max [45], are shown in Table 1.The compaction parameters were established using standard Proctor energy of 0.6 J/cm 3 .To prepare the material for further tests, an amount of water corresponding to w opt was admixed to the dry loess-sand mixture and left for an overnight rest in an airtight container.Two specimen sizes were used: 70 mm/140 mm and 38 mm/76 mm (diameter/height).Due to equipment limitations, slightly different procedures were applied to prepare specimens.The larger specimens were compacted in 3 layers with the standard Proctor energy directly in a split mold of appropriate dimensions.The smaller ones were cut out from the larger specimens by means of a cutting ring.Despite the different preparation procedures, the bulk densities of the specimens, irrespective of the size, differed by not more than 1%. the sand, loess, and the sand-loess mixture (after the removal of calcium carbonate; ISO 17892-4 [41]) are shown in Figure 2. The physical parameters of the mixture, namely sand (Sa), silt (Si), and clay (Cl) contents, mean grain size d50, coefficients of uniformity CU and curvature CC, plastic limit wP, liquid limit wL [42], plasticity index IP = wL − wP, colloidal activity A, loss on ignition LOI at 600 °C [43], particle density ρs [44], optimum moisture content wopt, and maximum dry density ρd.max [45], are shown in Table 1.The compaction parameters were established using standard Proctor energy of 0.6 J/cm 3 .To prepare the material for further tests, an amount of water corresponding to wopt was admixed to the dry loess-sand mixture and left for an overnight rest in an airtight container.Two specimen sizes were used: 70 mm/140 mm and 38 mm/76 mm (diameter/height).Due to equipment limitations, slightly different procedures were applied to prepare specimens.The larger specimens were compacted in 3 layers with the standard Proctor energy directly in a split mold of appropriate dimensions.The smaller ones were cut out from the larger specimens by means of a cutting ring.Despite the different preparation procedures, the bulk densities of the specimens, irrespective of the size, differed by not more than 1%.

Methodology
The prepared loess-sand specimens were subjected to unconfined compression tests, conventional triaxial compression tests, and bender-extender element tests in order to determine the shear strength and stiffness parameters of the material under various stress-strain and drainage conditions.A summary of all specimens tested is shown in Tables 2 and 3 with the initial values of the bulk density ρ0 and the void ratio e0.The following codes are used for the particular specimens: type of test/diameter for unconfined compression tests, e.g., UC/70, and type of test/diameter/confining stress for triaxial compression tests, e.g., CIU/38/100, CID/70/200.The triaxial tests conducted with the vertical stress paths are indicated with a letter V at the end, e.g., CID/70/200/V.
The unconfined compression tests were carried out, according to ISO 17892-7 [46] on six UC/38 and six UC/70 specimens.The specimens were unsaturated, with the average degree of saturation equal to 69%.All were compressed with an axial strain velocity of 25%/hour until a considerable drop in the vertical stress was observed.This velocity allowed for more than 10 readings to be taken before failure and provided the required testing time of between 2 and 15 min.The height change in the specimens was measured with an external LDT sensor.Additionally, LVDT on-specimen transducers measured the

Methodology
The prepared loess-sand specimens were subjected to unconfined compression tests, conventional triaxial compression tests, and bender-extender element tests in order to determine the shear strength and stiffness parameters of the material under various stressstrain and drainage conditions.A summary of all specimens tested is shown in Tables 2  and 3 with the initial values of the bulk density ρ 0 and the void ratio e 0 .The following codes are used for the particular specimens: type of test/diameter for unconfined compression tests, e.g., UC/70, and type of test/diameter/confining stress for triaxial compression tests, e.g., CIU/38/100, CID/70/200.The triaxial tests conducted with the vertical stress paths are indicated with a letter V at the end, e.g., CID/70/200/V.The unconfined compression tests were carried out, according to ISO 17892-7 [46] on six UC/38 and six UC/70 specimens.The specimens were unsaturated, with the average degree of saturation equal to 69%.All were compressed with an axial strain velocity of 25%/hour until a considerable drop in the vertical stress was observed.This velocity allowed for more than 10 readings to be taken before failure and provided the required testing time of between 2 and 15 min.The height change in the specimens was measured with an external LDT sensor.Additionally, LVDT on-specimen transducers measured the vertical strain and the change in the specimen diameter.The unconfined compression strength, the average E 50 modulus, and the Poisson's ratio ν were determined.
The conventional triaxial compression tests were performed according to ISO 17892-9 [47].The specimens were saturated at a back pressure of 450 kPa, achieving a Skempton's B value above 0.95.They were isotropically consolidated at the selected effective confining pressure σ' 3 (consolidation pressure σ' C ) in the range of 50-400 kPa.Next, the specimens were sheared, maintaining the cell pressure constant.Two specimens were subjected to a vertical stress path (constant p' = 200 or 350 kPa).The vertical displacement rate of the load frame was equal to 0.64%/h and 2.50%/h under the drained conditions (CID) and undrained (CIU) conditions, respectively.The shearing stage was continued until at least 15% of axial strain was achieved or until a clear shear failure surface was identified.The vertical deformation of the specimens was measured using external vertical displacement transducers.Local LVDT or Hall-effect sensors (c.f.Table 3) were additionally used on selected specimens.Unless stated otherwise, the strain calculations were based on the readings of the local displacement sensors, within their range.For larger strains and in the tests where the local sensors were not available, the readings of the external vertical displacement sensor and the volume change measurements were applied.The PE element tests (PETs) were performed on saturated and consolidated specimens 70 mm in diameter that were later subjected to drained shearing.Details of the PET methodology applied in this study can be found in [21].The readings were taken at σ' C (equal to the mean effective stress p' 0 ) within the range of 10-400 kPa.The triaxial compression and PETs were used to estimate the soil shear strength, as well as the values of the tangent and secant E and G moduli and their degradation with strain in drained and undrained conditions.The Poisson's ratio ν was also assessed.Various approaches were applied, as presented in Section 2.

Strength and Stiffness in Unconfined Compression
The average values of unconfined compression strength (UCS), secant Young's modulus (E 50 ), and Poisson's ratio (ν 50 ) (at the same strain for which E 50 was estimated) obtained in UC tests, together with the standard deviations (SD), are shown in Table 4.The values were calculated separately, based on the measurements from local (loc) and external (ex) displacement sensors.In Figure 3, showing the stress-strain curves for all the specimens tested in unconfined compression, quite a large scatter can be observed, which is common for this type of test [22].It is visually larger in UC/38 specimens (applying the same scale) and may result from the fact that smaller specimens are naturally less uniform, which is confirmed by larger SD values for the density and void ratio in Table 2.The differences in UCS obtained using loc and ex measurements for one specimen size result from the different area measurements and are negligible (<5%).The UCS values obtained are smaller than those reported for compacted natural loess (0.26-2.1 MPa) by other researchers [19,22,23].It shall be noted, however, that the mentioned results referred to the loess from Chinese Loess Plateau, which has a much higher natural calcium carbonate content, responsible for an increase in strength.In addition, the authors of [22] reported a large influence of moisture content on UCS values, and, therefore, the data cannot be directly compared.In this study, the effect of moisture content was not analyzed.However, the main factor responsible for the smaller UCS values is definitely the increased content of sand, resulting from the larger size of the voids and the smaller interparticle suction.The average UCS obtained for the UC/38 specimens is only 13% or 17% greater than for the UC/70 specimens, but it was achieved at much larger vertical strain (ex: 5.5-7.5% vs. 1.5-2.5% or loc: 5.0-7.5% vs. 1.0-1.1%),eventually giving much smaller Young's modulus values.This may result from both different sample sizes and different preparation procedures and requires further investigation.The average E 50 moduli of the UC/38 specimens seem to not depend much on the method of measurement (loc or ex), while in the UC/70 specimens the stiffness measured with the external LDT sensor is much smaller than the one based on the onspecimen transducers.As the tests of different specimen sizes were performed also in different apparatuses, this may be connected with the different sensitivity of the sensors used or the typical problems encountered in conventional deformation measurements, such as seating or bedding errors [48,49].
showed a large scatter; thus, only the four most representative specimens have been shown.Regardless of the size of the specimen, the Poisson's ratio increases gradually with SF.A similar phenomenon was reported by Bowers [24] in triaxial tests on unsaturated loess.At SF = 0.5 (denoted as ν50), the Poisson's ratio is equal to approximately 0.14 (for both sizes), which is consistent with the results presented in [50] and close to 0.1, which is a value typically assumed for unsaturated soils [51].

Shear Strength in Triaxial Compression
The stress paths p'-q and the shearing characteristics ε1-q obtained in the triaxial compression tests are shown in Figures 5 and 6.The circular and square markers indicate the points representing the criteria (q/p')max and (q)max, respectively.For clarity, they were presented only for CIU/38 and CID/38 specimens.It is clear that the behavior of the soil depends strongly on the drainage conditions, which is mainly due to the different shapes of the stress paths.In the drained tests, the peak shear strength (q)max is much smaller than The Poisson's ratio ν was calculated based on the local strain measurements.Figure 4 presents the variation of ν with stress ratio SF = σ v /σ v.max .The results for UC/70 specimens showed a large scatter; thus, only the four most representative specimens have been shown.Regardless of the size of the specimen, the Poisson's ratio increases gradually with SF.A similar phenomenon was reported by Bowers [24] in triaxial tests on unsaturated loess.At SF = 0.5 (denoted as ν 50 ), the Poisson's ratio is equal to approximately 0.14 (for both sizes), which is consistent with the results presented in [50] and close to 0.1, which is a value typically assumed for unsaturated soils [51].
different specimen sizes were performed also in different apparatuses, this may be connected with the different sensitivity of the sensors used or the typical problems encountered in conventional deformation measurements, such as seating or bedding errors [48,49].The Poisson's ratio  was calculated based on the local strain measurements.Figure 4 presents the variation of  with stress ratio SF = σv/σv.max.The results for UC/70 specimens showed a large scatter; thus, only the four most representative specimens have been shown.Regardless of the size of the specimen, the Poisson's ratio increases gradually with SF.A similar phenomenon was reported by Bowers [24] in triaxial tests on unsaturated loess.At SF = 0.5 (denoted as ν50), the Poisson's ratio is equal to approximately 0.14 (for both sizes), which is consistent with the results presented in [50] and close to 0.1, which is a value typically assumed for unsaturated soils [51].

Shear Strength in Triaxial Compression
The stress paths p'-q and the shearing characteristics ε1-q obtained in the triaxial compression tests are shown in Figures 5 and 6.The circular and square markers indicate the points representing the criteria (q/p')max and (q)max, respectively.For clarity, they were presented only for CIU/38 and CID/38 specimens.It is clear that the behavior of the soil depends strongly on the drainage conditions, which is mainly due to the different shapes of the stress paths.In the drained tests, the peak shear strength (q)max is much smaller than

Shear Strength in Triaxial Compression
The stress paths p'-q and the shearing characteristics ε 1 -q obtained in the triaxial compression tests are shown in Figures 5 and 6.The circular and square markers indicate the points representing the criteria (q/p') max and (q) max , respectively.For clarity, they were presented only for CIU/38 and CID/38 specimens.It is clear that the behavior of the soil depends strongly on the drainage conditions, which is mainly due to the different shapes of the stress paths.In the drained tests, the peak shear strength (q) max is much smaller than in the undrained tests.It is also achieved at smaller axial strains: about 1.1-4.0%and 4.5-6.9% in the CID/70 and CID/38 specimens, compared to >15% and 7.7-16.0%in the CIU/70 and CIU/38 specimens, respectively.For drained tests, the peak corresponds to the failure condition (q/p') max , while in undrained tests, the (q/p') max is achieved at significantly lower strains-approximately 0.8-1.5% and 1.7-5.7% in the CIU/70 and CIU/38 tests, respectively.The strain at which the failure boundary is reached is, therefore, affected not only by the type of drainage but also by the size of the specimen.Furthermore, the stress-strain curves show more distinct peaks in drained conditions and in the smaller specimens.The influence of the size of the specimen on the softening rate was also observed by Hu et al. [52] and may be related to the fact that in the larger specimens more than one distinct failure surface was observed, causing a different stress redistribution.The size of the specimen also played a role in the evolution of excess pore pressure u in the undrained tests, contrary to the observations of Omar and Sadrekarimi [53] in their tests on sand.For both specimen sizes, the excess pore pressure initially increased with axial strain, after which it started to decrease, reaching negative values (suction).For large specimens, the maximum u values were obtained at strains of approximately 0.4-0.8%,that is, before (q/p') max was reached, and negative u values were obtained at strains of approximately 1.2-3.8%,that is, after (q/p') max and long before (q) max were achieved.For small specimens, usually both the maximum and negative values appeared before reaching the (q/p') max , that is, at about 0.3-1.5% and 1.3-6.2%,respectively.For drained tests, after an initial decrease in volume with the peak at approximately 0.2-1.6%(CID/70) or 2.0-3.5% (CID/38) of axial strains, the attainment of negative volumetric strain values occurred approximately when the failure condition was reached ((q) max ~(q/p') max ), irrespective of the specimen size.It is worth mentioning that the strain evolution of the CID/70/350/V and CIU/70/400 specimens, as well as the CID/70/200/V and CIU/70/200 specimens, with very similar stress paths, was practically the same until the stress path reached the failure boundary; then, under drained conditions, the specimens began to yield, while in the undrained tests, as the stress path 'slid' along the failure boundary, a constant increase in shear strength was observed.This strain hardening phenomenon is typical for strongly dilatant well-compacted loess, cf.e.g., [12,15].     5.They were calculated based on two selected failure criteria: (q/p') max and (q) max .The steady or critical state, in which the mass is continuously deforming at constant volume, constant normal effective stress, constant shear stress, and constant velocity [54], was not observed until 15% axial strain; therefore, this criterion was not considered.It should be noted that within the stress range analyzed, the FBs considered were rectilinear, with a coefficient of determination of R 2 ≥ 0.999.The selection of the failure criterion only had an influence on the results of the CIU tests, as, in the drained tests, both criteria occurred at almost the same strains.With the (q) max considered as the shear strength, slightly smaller φ' values were obtained at (forced) zero cohesion.As only very slight strength degradation and pore pressure change were observed in these tests after (q) max was reached, these results are considered to be close to the steady or critical state.There is no clear dependency between the angle of internal friction and the size of the specimen, but it appears that, except for the undrained test at the (q/p') max criterion, φ' becomes slightly smaller (by 1.4-2.4• ) in the larger specimens, which is similar to the trends observed in sands [52,53].Cohesion also tends to decrease in the larger specimens.The angle of internal friction values are high, when compared to the results obtained by other researchers for natural compacted loess, which have been summarized in Table 6.This is mainly due to the lower clay content and the higher sand content in the loess-sand mixture tested in this study, indicating the effectiveness of the improvement in terms of the angle of internal friction.The relatively lower cohesion can be explained by the smaller amount of clay fraction (and so clayey minerals), which is believed to coat the individual silt grains and to be the main cementing agent in the natural loess [55], next to calcite.It should be noted that the shear strength parameters obtained in this investigation correspond well to the values obtained by Xu et al. [15] for the loess with similar physical properties.The results of the initial shear moduli G 0 , determined in the PETs, in relation to the mean effective stress p' 0 , are given in Figure 7a (after [21]).The G 0 values were compared with the results of three other studies on similar soils: (a) by Wang et al. [13] on normally consolidated and over-consolidated compacted loess (Sa = 5%, Si = 86%, Cl = 9%; initial bulk density ρ 0 = 1.7 g/cm 3 ), (b) by Song et al. [20] on undisturbed Lanzhou loess wetted to w ≈ 16% (Sa = 5%, Si = 85%, Cl = 10%, ρ 0 = 1.45-1.75g/cm 3 ), and (c) an empirical equation proposed by Hassanipour et al. [58] for well-compacted sand-clay mixtures with Sa ≤ 60%.The G 0 values of the loess-sand mixture in our research were almost twice as large as of the loess with no admixtures reported in the literature, indicating the strong positive influence of compaction and increased sand content.They were, however, lower than in the sand-clay mixture.The E0 values calculated with Equation ( 8) are shown in Figure 7b (after [21]).They can be roughly estimated as being equal to 2.97 G0.The dependence between the Young's modulus and p'0 may be described with the following Equation (11):  The dependence between the initial stiffness modulus of the loess-sand mixture and p' 0 can be expressed with the use of a simple power law widely accepted in soil dynamics, resulting in the following Equation (10): (10) in which p a is the atmospheric pressure equal to 101 kPa.According to Equation ( 8), to evaluate the initial elastic modulus E 0 , it is necessary to know the value of the Poisson's ratio ν 0 .The values obtained from the PETs, based on Equation ( 9), are presented in Figure 8, showing no significant dependence on the mean effective stress.The average v value is equal to 0.485 ± 0.03, indicating nearly undrained conditions during PE testing.In fully saturated specimens, the first observed P-wave belongs to water travelling with the velocity V P ≈ 1485 m/s.The measured P-wave velocities in the loess-sand mixtures increased from ca. 950 m/s at p' 0 = 10 kPa to ca. 1980 m/s at p' 0 = 400 kPa.They reached 1485 m/s at p' 0 = 100 kPa.The V P values lower than 1485 m/s at the smallest confining stresses may result from a weaker contact between the piezoceramic element and soil.The values of the Poisson's ratio obtained in the PE tests were much higher than in the UC tests, which can be explained with the much higher degree of saturation S. The increase in the value of the Poisson's ratio with S is known from the literature [50,59].
piezoceramic element and soil.The values of the Poisson's ratio obtained in the PE tests were much higher than in the UC tests, which can be explained with the much higher degree of saturation S. The increase in the value of the Poisson's ratio with S is known from the literature [50,59].
The E0 values calculated with Equation ( 8) are shown in Figure 7b (after [21]).They can be roughly estimated as being equal to 2.97 G0.The dependence between the Young's modulus and p'0 may be described with the following Equation ( 11 9)) with the mean effective stress p'0.

Shear Modulus G
It is well known that the shear stiffness of soil decreases with shear strain γ.Such characteristics were also obtained for the loess-sand mixture tested.The degradation curves of the secant and tangent shear moduli, Gs and Gt, are presented in Figures 9 and  10.For clarity, only the CIU/38 and CID/38 specimens are included in these graphs.To eliminate the influence of the mean effective stress, the results have been normalized with G0, as calculated from Equation (10).As can be seen, the applied procedure was very effective in reducing the scatter of data, at least for the secant modulus.The E 0 values calculated with Equation ( 8) are shown in Figure 7b (after [21]).They can be roughly estimated as being equal to 2.97 G 0 .The dependence between the Young's modulus and p' 0 may be described with the following Equation ( 11):

Shear Modulus G
It is well known that the shear stiffness of soil decreases with shear strain γ.Such characteristics were also obtained for the loess-sand mixture tested.The degradation curves of the secant and tangent shear moduli, G s and G t , are presented in Figures 9 and 10.For clarity, only the CIU/38 and CID/38 specimens are included in these graphs.To eliminate the influence of the mean effective stress, the results have been normalized with G 0 , as calculated from Equation (10).As can be seen, the applied procedure was very effective in reducing the scatter of data, at least for the secant modulus.There are many theoretical and empirical approaches to analytically describe the non-linear stress-strain behavior of soils [60].The G/G0 degradation curve is most often modelled with a hyperbole, as first proposed by Hardin and Drnevich [61].One of the most commonly used modifications of that proposal has been given by Darendeli [62], as follows: where γref is the reference shear strain and α is the curvature coefficient.Empirical equations and some further extensions based on this formula can be found, e.g., in [63- There are many theoretical and empirical approaches to analytically describe the non-linear stress-strain behavior of soils [60].The G/G0 degradation curve is most often modelled with a hyperbole, as first proposed by Hardin and Drnevich [61].One of the most commonly used modifications of that proposal has been given by Darendeli [62], as follows: where γref is the reference shear strain and α is the curvature coefficient.Empirical equations and some further extensions based on this formula can be found, e.g., in [63- There are many theoretical and empirical approaches to analytically describe the non-linear stress-strain behavior of soils [60].The G/G 0 degradation curve is most often modelled with a hyperbole, as first proposed by Hardin and Drnevich [61].One of the most commonly used modifications of that proposal has been given by Darendeli [62], as follows: where γ ref is the reference shear strain and α is the curvature coefficient.Empirical equations and some further extensions based on this formula can be found, e.g., in [63][64][65].For fine soils, researchers often include the plasticity index I P as a variable in the equation.This type of G s /G 0 model (Equation ( 12)) was verified by Vardanega and Bolton [66] in a large set of the literature results for clays and silts.In static applications, they achieved the best fit at α = 0.736 and γ ref = 0.0022•I P (where I P is expressed numerically, not as a percentage).
In this study, the same approach with I P = 0.0564 was applied, and the results are presented in Figures 9 and 10 as black dashed lines 'V&B (2013) IP'.In terms of shape, the empirical equation shows quite a good agreement with the results of both drained and undrained secant degradation curves, although quantitatively it shall rather be treated as a lower boundary of the G s /G 0 values.This non-compliance may result from very low colloidal activity (typical for loess).For fine textbook soils, the colloidal activity is around 1, meaning that the plasticity index is close to their clay content.For loess, I P is significantly higher than the content of the clay fraction and, thus, the plasticity index may be less descriptive of the soil type.In the analyzed case, the use of clay content instead of the plasticity index in the γ ref calculation gave a better fit for the degradation curves, which is shown with the green dashed lines 'V&B (2013) Cl'.The results obtained in this study were also compared with the model for sands presented by Oztoprak and Bolton [64], in which an additional parameter γ e was introduced, shifting the degradation curve rightwards in order to describe an extended zone of elastic behavior, as follows: If γ e = 0 is assumed, then Equation ( 13) reduces to Equation (12).Equation ( 13), with the parameters representing an average sand (γ ref = 0.044%, γ e = 0.0007%, and α = 0.88) was used to plot the dotted lines 'O&B (2013)' in Figures 9 and 10.The model of Oztoprak and Bolton [64] may be considered an upper boundary for the secant moduli results obtained.
Other empirical models available in the literature gave much weaker fitting to the laboratory data, and so are not discussed here.The hyperbolic models describe the course of degradation curve quite well in the small strain region.However, they are not sufficient to fully characterize the complex behavior of soils in the full range of strains.The discrepancy between the theoretical curves and the readings increases at shear strains greater than 1%, as presented in Figure 11.The influence of the drainage conditions becomes clear as well.In undrained tests, between shear strains of approximately 3-11% (which correspond to approximately 2-8% of the axial strain), a flattening of the degradation curves is visible.In drained tests, this effect is not apparent, and the stiffness degradation continues with increasing strain.Nevertheless, irrespective of the drainage conditions, the G s /G 0 values at large shear strains are much higher than the theoretical predictions.However, it should be noted that the theoretical models invoked were based on different types of tests (resonant column, triaxial test, torsional shear and direct shear) and procedures with different strain range and the influence of drainage conditions was not explored in detail.The shear strain range at which the curve levelling is observed in CIU tests is not directly related to the occurrence of the (q) max or (q/p) max failure conditions.However, it seems to be connected with the cross-coupling between volumetric and distortional effects, cf.[30].In the CID tests, the specimen initially becomes stiffer as it contracts and softens at larger strains, i.e., as it dilates, while in the undrained test, no volume changes are allowed.Due to the fact that researchers rarely perform simultaneous tests with and without drainage, it is difficult to support these considerations with the results of other studies.The explanation of this phenomenon requires further investigation.The tangent moduli obtained experimentally are smaller than the secant moduli (cf.Figures 9 and 10), as expected [67].Up to γ~1%, the degradation curves Gt/G0 are located below the Gs/G0 curves and also below the empirical models that were examined.When strain softening is observed, the tangent moduli become negative, although the absolute values of Gt/G0 are very small.Like for the secant moduli, at shear strains greater than about 1%, the behavior depends on the drainage conditions.The difference between the values of the tangent and secant moduli obviously results from the different definitions and calculation methods.The tangent modulus is very sensitive to the smallest changes in the readings-this is why a larger scatter of values occurs for Gt/G0, when compared to Gs/G0.Additionally, the value of the tangent modulus depends on the range of data taken into consideration (the i−1 and i+1 readings) and, thus, on the frequency of the readings.There is no unambiguous relationship between the tangent and secant moduli, as presented in Figure 12.The mentioned problems are the most probable reasons why in the available literature there are much fewer data for secant moduli than for tangent moduli, and no empirical models are available to describe this degradation.The tangent moduli obtained experimentally are smaller than the secant moduli (cf.Figures 9 and 10), as expected [67].Up to γ~1%, the degradation curves G t /G 0 are located below the G s /G 0 curves and also below the empirical models that were examined.When strain softening is observed, the tangent moduli become negative, although the absolute values of G t /G 0 are very small.Like for the secant moduli, at shear strains greater than about 1%, the behavior depends on the drainage conditions.The difference between the values of the tangent and secant moduli obviously results from the different definitions and calculation methods.The tangent modulus is very sensitive to the smallest changes in the readings-this is why a larger scatter of values occurs for G t /G 0 , when compared to G s /G 0 .Additionally, the value of the tangent modulus depends on the range of data taken into consideration (the i−1 and i+1 readings) and, thus, on the frequency of the readings.There is no unambiguous relationship between the tangent and secant moduli, as presented in Figure 12.The mentioned problems are the most probable reasons why in the available literature there are much fewer data for secant moduli than for tangent moduli, and no empirical models are available to describe this degradation.The tangent moduli obtained experimentally are smaller than the secant moduli (cf.Figures 9 and 10), as expected [67].Up to γ~1%, the degradation curves Gt/G0 are located below the Gs/G0 curves and also below the empirical models that were examined.When strain softening is observed, the tangent moduli become negative, although the absolute values of Gt/G0 are very small.Like for the secant moduli, at shear strains greater than about 1%, the behavior depends on the drainage conditions.The difference between the values of the tangent and secant moduli obviously results from the different definitions and calculation methods.The tangent modulus is very sensitive to the smallest changes in the readings-this is why a larger scatter of values occurs for Gt/G0, when compared to Gs/G0.Additionally, the value of the tangent modulus depends on the range of data taken into consideration (the i−1 and i+1 readings) and, thus, on the frequency of the readings.There is no unambiguous relationship between the tangent and secant moduli, as presented in Figure 12.The mentioned problems are the most probable reasons why in the available literature there are much fewer data for secant moduli than for tangent moduli, and no empirical models are available to describe this degradation.It was found that in the drained tests, at the shear strains of up to 3-4%, the G s /G 0 values were generally larger (slower stiffness degradation) in the 70 mm diameter specimens than in the smaller ones, which is visible in Figure 13b.This is consistent with the observations from the UC tests, cf. the E 50 values in Table 4.It can be noticed that the curves CID/70/100 and CID/70/200 are of different shapes than the curves CID/70/50 and CID/70/350, which results from the fact that in these tests the readings of the local displacement sensors were not available.In the CID/38 tests, the differences between the external and local strain readings were much smaller, and no difference in the shape of the curve was found.In the undrained tests (Figure 13a), up to 10% shear strain, the moduli of the larger specimens were also higher than the ones of smaller specimen.No flattening effect was observed in the 70 mm-diameter specimens.However, the CIU/70/50 specimen was exceptional.As the stiffness calculated on the basis of external sensors is usually smaller than that obtained using local ones, the shear moduli in the CIU/70 series would probably be even higher, but the trend is comparable.The number of specimens examined is too small to draw general conclusions, but the data presented indicate that the specimen size can have a significant effect on the values of the shear moduli and their strain degradation.Furthermore, evaluation of the degradation curves based on the external displacement measurements shall be avoided due to the possible seating and bedding errors, as mentioned earlier.It was found that in the drained tests, at the shear strains of up to 3-4%, the Gs/G0 values were generally larger (slower stiffness degradation) in the 70 mm diameter specimens than in the smaller ones, which is visible in Figure 13b.This is consistent with the observations from the UC tests, cf. the E50 values in Table 4.It can be noticed that the curves CID/70/100 and CID/70/200 are of different shapes than the curves CID/70/50 and CID/70/350, which results from the fact that in these tests the readings of the local displacement sensors were not available.In the CID/38 tests, the differences between the external and local strain readings were much smaller, and no difference in the shape of the curve was found.In the undrained tests (Figure 13a), up to 10% shear strain, the moduli of the larger specimens were also higher than the ones of smaller specimen.No flattening effect was observed in the 70 mm-diameter specimens.However, the CIU/70/50 specimen was exceptional.As the stiffness calculated on the basis of external sensors is usually smaller than that obtained using local ones, the shear moduli in the CIU/70 series would probably be even higher, but the trend is comparable.The number of specimens examined is too small to draw general conclusions, but the data presented indicate that the specimen size can have a significant effect on the values of the shear moduli and their strain degradation.Furthermore, evaluation of the degradation curves based on the external displacement measurements shall be avoided due to the possible seating and bedding errors, as mentioned earlier.A significant difference was observed between the stiffness degradation curves for the specimens sheared along conventional drained (σ'3 = σ'C = const) and vertical (p' = σ'C = const) stress paths.The results are presented in Figure 14 with circles and triangles, respectively.It is worth noting that the G0 values used for normalization (207.9MPa for 'C = 200 kPa and 290.1 MPa for 'C = 350 kPa) were obtained at the isotropic state of stress, corresponding to the beginning of the paths.Both the secant and the tangent moduli obtained for vertical stress paths are smaller than the moduli from conventional compression tests.This can be explained by the fact that at the same q the ratio q/qmax is greater in the tests with the vertical stress path, due to the much smaller values of qmax achieved.Therefore, a greater shear strain develops, and the stiffness modulus is smaller.A significant difference was observed between the stiffness degradation curves for the specimens sheared along conventional drained (σ' 3 = σ' C = const) and vertical (p' = σ' C = const) stress paths.The results are presented in Figure 14 with circles and triangles, respectively.It is worth noting that the G 0 values used for normalization (207.9MPa for σ' C = 200 kPa and 290.1 MPa for σ' C = 350 kPa) were obtained at the isotropic state of stress, corresponding to the beginning of the paths.Both the secant and the tangent moduli obtained for vertical stress paths are smaller than the moduli from conventional compression tests.This can be explained by the fact that at the same q the ratio q/q max is greater in the tests with the vertical stress path, due to the much smaller values of q max achieved.Therefore, a greater shear strain develops, and the stiffness modulus is smaller.The Poisson's ratio ν is often assumed to be constant.It is obviously true in undrained triaxial tests, but not in the tests with drainage allowed.This can be seen in Figure 15, where the ν values of CID/38 specimens, calculated from Equation (3), are plotted in relation to the stress ratio SF = q/(q)max.The Poisson's ratios are initially close to the ν0 values from the BE tests and gradually decrease to 0.35-0.40 as the stress intensity increases.However, starting from SF = 0.5, ν increases again, reaching values between 0.5-0.6 after the peak shear stress.A Poisson's ratio greater than 0.5 is treated as a nonphysical value, which means that the volume of the specimen increases during shearing.Actually, this is true for strongly dilatant soils, such as the well compacted loess-sand mixture under investigation.In general, it is assumed that the degradation behavior of the Young's modulus is the same as that of the shear modulus [68].This holds true only if the Poisson's ratio does not The Poisson's ratio ν is often assumed to be constant.It is obviously true in undrained triaxial tests, but not in the tests with drainage allowed.This can be seen in Figure 15, where the ν values of CID/38 specimens, calculated from Equation (3), are plotted in relation to the stress ratio SF = q/(q) max .The Poisson's ratios are initially close to the ν 0 values from the BE tests and gradually decrease to 0.35-0.40 as the stress intensity increases.However, starting from SF = 0.5, ν increases again, reaching values between 0.5-0.6 after the peak shear stress.A Poisson's ratio greater than 0.5 is treated as a nonphysical value, which means that the volume of the specimen increases during shearing.Actually, this is true for strongly dilatant soils, such as the well compacted loess-sand mixture under investigation.The Poisson's ratio ν is often assumed to be constant.It is obviously true in undrained triaxial tests, but not in the tests with drainage allowed.This can be seen in Figure 15, where the ν values of CID/38 specimens, calculated from Equation (3), are plotted in relation to the stress ratio SF = q/(q)max.The Poisson's ratios are initially close to the ν0 values from the BE tests and gradually decrease to 0.35-0.40 as the stress intensity increases.However, starting from SF = 0.5, ν increases again, reaching values between 0.5-0.6 after the peak shear stress.A Poisson's ratio greater than 0.5 is treated as a nonphysical value, which means that the volume of the specimen increases during shearing.Actually, this is true for strongly dilatant soils, such as the well compacted loess-sand mixture under investigation.In general, it is assumed that the degradation behavior of the Young's modulus is the same as that of the shear modulus [68].This holds true only if the Poisson's ratio does not In general, it is assumed that the degradation behavior of the Young's modulus is the same as that of the shear modulus [68].This holds true only if the Poisson's ratio does not vary significantly during the loading process.This is not the case for the drained tests, as shown above.In the literature, Young's modulus E is more often normalized with p' 0 [67,69] than with E 0 [34,70] or plotted against the stress ratio SF [68,71].This is probably due to the difficulty in determining the initial modulus E 0 in triaxial tests, resulting from the fact that the measurement of the P-wave velocity is much less available and is more troublesome than the S-wave.In this study, the values of the Young's modulus were normalized with E 0 , calculated from Equation ( 8), to maintain a methodology similar to that used for the shear modulus.The degradation curves of the secant and tangent elastic moduli with axial strain ε 1 for 38 mm diameter specimens are presented in Figures 16 and 17, respectively.Some theoretical relations describing the E/E 0 degradation curves can be found in [34,70,72].However, there are practically no results for loess available in the literature to compare.Like for G, the tangent Young's modulus is smaller than the shear modulus.The effect of the specimen size was also observed, as well as the influence of drainage conditions on the degradation curve at strains greater than 1%.
Materials 2024, 17, x FOR PEER REVIEW 18 of 23 vary significantly during the loading process.This is not the case for the drained tests, as shown above.In the literature, Young's modulus E is more often normalized with p'0 [67,69] than with E0 [34,70] or plotted against the stress ratio SF [68,71].This is probably due to the difficulty in determining the initial modulus E0 in triaxial tests, resulting from the fact that the measurement of the P-wave velocity is much less available and is more troublesome than the S-wave.In this study, the values of the Young's modulus were normalized with E0, calculated from Equation ( 8), to maintain a methodology similar to that used for the shear modulus.The degradation curves of the secant and tangent elastic moduli with axial strain ε1 for 38 mm diameter specimens are presented in Figure 16 and Figure 17, respectively.Some theoretical relations describing the E/E0 degradation curves can be found in [34,70,72].However, there are practically no results for loess available in the literature to compare.Like for G, the tangent Young's modulus is smaller than the shear modulus.The effect of the specimen size was also observed, as well as the influence of drainage conditions on the degradation curve at strains greater than 1%.vary significantly during the loading process.This is not the case for the drained tests, as shown above.In the literature, Young's modulus E is more often normalized with p'0 [67,69] than with E0 [34,70] or plotted against the stress ratio SF [68,71].This is probably due to the difficulty in determining the initial modulus E0 in triaxial tests, resulting from the fact that the measurement of the P-wave velocity is much less available and is more troublesome than the S-wave.In this study, the values of the Young's modulus were normalized with E0, calculated from Equation ( 8), to maintain a methodology similar to that used for the shear modulus.The degradation curves of the secant and tangent elastic moduli with axial strain ε1 for 38 mm diameter specimens are presented in Figure 16 and Figure 17, respectively.Some theoretical relations describing the E/E0 degradation curves can be found in [34,70,72].However, there are practically no results for loess available in the literature to compare.Like for G, the tangent Young's modulus is smaller than the shear modulus.The effect of the specimen size was also observed, as well as the influence of drainage conditions on the degradation curve at strains greater than 1%.As mentioned at the beginning, E 50 is the specific secant modulus that is often required in practical applications and in numerical modelling.The E 50 values, together with the Poisson's ratios v 50 , obtained at the same strain as E 50 , are shown in Table 7.The moduli and Poisson's ratios obtained in drained and undrained tests are denoted as E' 50 , v' 50 , E u50 , and v u50 , respectively.The relation between E 50 and consolidation pressure σ' C is also presented in Figure 18.For comparison, the values obtained in UC tests are shown at σ' C = 0 kPa, although the true effective confining stress (suction) is unknown.The relation can be described using a power function, with a reasonably good fit.It is evident that the larger specimens showed higher E 50 values, both in tests with and without drainage, which is consistent with the results of the UC tests.Greater E 50 values were generally obtained in the drained tests, which is related to the lower axial strains at which the stress intensity was equal to 50% of (q) max ; this is visible in Figure 17, where E 50 values have been marked with 'x' on the CID/38 and CIU/38 degradation curves.Only for the 38 mm in diameter specimens, at σ' C ≤ 150 kPa, were the E 50 values slightly higher in the undrained than in the drained conditions.The results obtained in CIU tests are similar to those reported by Zarei et al. [73] for specimens of a compacted loess, 50 mm in diameter and 100 mm height, with very similar density.Contrary to the E 50 modulus, which for the loess-sand mixture depends on both specimen size and consolidation pressure, the Poisson's ratio ν 50 can be assumed as constant and equal to 0.5 and 0.38 (±0.03) in the undrained and drained conditions, respectively.As mentioned at the beginning, E50 is the specific secant modulus that is often required in practical applications and in numerical modelling.The E50 values, together with the Poisson's ratios v50, obtained at the same strain as E50, are shown in Table 7.The moduli and Poisson's ratios obtained in drained and undrained tests are denoted as E'50, v '50, Eu50, and vu50, respectively.The relation between E50 and consolidation pressure 'C is also presented in Figure 18.For comparison, the values obtained in UC tests are shown at 'C = 0 kPa, although the true effective confining stress (suction) is unknown.The relation can be described using a power function, with a reasonably good fit.It is evident that the larger specimens showed higher E50 values, both in tests with and without drainage, which is consistent with the results of the UC tests.Greater E50 values were generally obtained in the drained tests, which is related to the lower axial strains at which the stress intensity was equal to 50% of (q)max; this is visible in Figure 17, where E50 values have been marked with 'x' on the CID/38 and CIU/38 degradation curves.Only for the 38 mm in diameter specimens, at σ'C ≤ 150 kPa, were the E50 values slightly higher in the undrained than in the drained conditions.The results obtained in CIU tests are similar to those reported by Zarei et al. [73] for specimens of a compacted loess, 50 mm in diameter and 100 mm height, with very similar density.Contrary to the E50 modulus, which for the loess-sand mixture depends on both specimen size and consolidation pressure, the Poisson's ratio ν50 can be assumed as constant and equal to 0.5 and 0.38 (±0.03) in the undrained and drained conditions, respectively.

Summary and Conclusions
The parameters commonly used to describe the shear strength and stiffness of soils vary depending on the testing conditions.In this study, a compacted loess with 20% sand admixture was investigated in monotonic triaxial and unconfined compression tests.Its shear strength was evaluated using different criteria.The values of shear and elastic moduli and the Poisson's ratio have been determined, together with their dependence on strain and stress range, draining conditions, and specimen size.Both tangent and secant values of the stiffness moduli were calculated.The main conclusions are as follows: • The UCS values of the loess-sand mixture were not much influenced by the size of the specimen.They were smaller than reported by other researchers for naturally compacted loess, which is due to the higher sand content and, possibly, the lower CaCO 3 content.• The choice of the shear strength criterion ((q) max or (q/p) max ), the specimen size and drainage conditions had little influence on the value of effective angle of internal friction φ', which was slightly larger than that of typical natural loess.On the other hand, the values of cohesion c' depended on both the criterion used and the size of the specimens.
• The initial moduli G 0 and E 0 increase with consolidation pressure, and the dependency can be described with a classical power function.The values obtained for the loesssand mixture are higher than those published for natural loess, due to the higher density and sand content, confirming the positive influence of compaction and sand admixture on loess stiffness.• Tangent stiffness moduli show higher variability, as they are very sensitive to any changes in the slope of the stress-strain curve.Their values are smaller than the secant moduli, and after achieving maximum stress intensity, they fall below zero.There is no obvious relationship between the secant and tangent moduli.• The larger specimens have usually been found to be stiffer.The shape of the stress path also influences the soil stiffness, which is related to the current stress level.The specimens sheared along vertical effective p'-q stress paths show faster degradation of stiffness with strain.• The clay content is more representative for loess than I P for the theoretical description of stiffness degradation.• The Poisson's ratio ν depends on the degree of saturation and is not constant during drained shearing.No dependency on the size of the specimen was identified.• The Young's modulus E shows behavior similar to that of the shear modulus.The E 50 values are greater in larger specimens, regardless of the confining stress.They are also greater in drained than in undrained triaxial tests, which is mainly due to the much smaller strain at which 50% ultimate stress intensity is achieved.Due to the high sensitivity of E 50 to the testing conditions, it is not recommended as a valuable parameter in numerical analyses.

Figure 3 .Figure 4 .
Figure 3. Shear characteristics for six UC/38 (a) and six UC/70 (b) specimens based on external sensors (black dashed lines) and local sensors (red continuous lines).

Figure 3 .
Figure 3. Shear characteristics for six UC/38 (a) and six UC/70 (b) specimens based on external sensors (black dashed lines) and local sensors (red continuous lines).

Figure 3 .Figure 4 .
Figure 3. Shear characteristics for six UC/38 (a) and six UC/70 (b) specimens based on external sensors (black dashed lines) and local sensors (red continuous lines).
Materials 2024, 17, x FOR PEER REVIEW 10 of 23 parameters obtained in this investigation correspond well to the values obtained by Xu et al. [15] for the loess with similar physical properties.

Figure 6 .
Figure 6.Stress paths (a) and shear characteristics (b) for CID/38 and CID/70 specimens.The effective strength parameters (cohesion c' and angle of internal friction φ') of the Mohr-Coulomb failure boundary (FB) are summarized in Table5.They were calculated based on two selected failure criteria: (q/p') max and (q) max .The steady or critical state, in which the mass is continuously deforming at constant volume, constant normal effective stress, constant shear stress, and constant velocity[54], was not observed until 15% axial strain; therefore, this criterion was not considered.It should be noted that within the stress range analyzed, the FBs considered were rectilinear, with a coefficient of determination of R 2 ≥ 0.999.The selection of the failure criterion only had an influence on the results of the CIU tests, as, in the drained tests, both criteria occurred at almost the same strains.With the (q) max considered as the shear strength, slightly smaller φ' values were obtained at (forced) zero cohesion.As only very slight strength degradation and pore pressure change were observed in these tests after (q) max was reached, these results are considered to be close to the steady or critical state.There is no clear dependency between the angle of internal friction and the size of the specimen, but it appears that, except for the undrained test at the (q/p') max criterion, φ' becomes slightly smaller (by 1.4-2.4• ) in the larger specimens, which is similar to the trends observed in sands[52,53].Cohesion also tends to decrease in the larger specimens.The angle of internal friction values are high, when compared to the results obtained by other researchers for natural compacted loess, which have been summarized in Table6.This is mainly due to the lower clay content and the higher sand content in the loess-sand mixture tested in this study, indicating the effectiveness of the improvement in terms of the angle of internal friction.The relatively lower cohesion can be explained by the smaller amount of clay fraction (and so clayey minerals), which is believed to coat the individual silt grains and to be the main cementing agent in the natural loess[55], next to calcite.It should be noted that the shear strength parameters obtained in this investigation correspond well to the values obtained by Xu et al.[15] for the loess with similar physical properties.

Materials 2024 ,
17, x FOR PEER REVIEW 12 of 23velocities in the loess-sand mixtures increased from ca. 950 m/s at p'0 = 10 kPa to ca. 1980 m/s at p'0 = 400 kPa.They reached 1485 m/s at p'0 = 100 kPa.The VP values lower than 1485 m/s at the smallest confining stresses may result from a weaker contact between the piezoceramic element and soil.The values of the Poisson's ratio obtained in the PE tests were much higher than in the UC tests, which can be explained with the much higher degree of saturation S. The increase in the value of the Poisson's ratio with S is known from the literature[50,59].

Figure 7 .
Figure 7. Influence of the mean effective stress p'0 on (a) the initial shear modulus G0 and (b) the initial Young's modulus E0 [9].Figure 7. Influence of the mean effective stress p' 0 on (a) the initial shear modulus G 0 and (b) the initial Young's modulus E 0 [9].

Figure 7 .
Figure 7. Influence of the mean effective stress p'0 on (a) the initial shear modulus G0 and (b) the initial Young's modulus E0 [9].Figure 7. Influence of the mean effective stress p' 0 on (a) the initial shear modulus G 0 and (b) the initial Young's modulus E 0 [9].

Figure 7 .
Figure 7. Influence of the mean effective stress p'0 on (a) the initial shear modulus G0 and (b) the initial Young's modulus E0 [9].

Figure 13 .
Figure 13.Comparison of the degradation curves Gs/G0 for specimens with different sizes (38 and 70 mm in diameter) in undrained CIU (a) and drained CID (b) triaxial tests.

Figure 13 .
Figure 13.Comparison of the degradation curves G s /G 0 for specimens with different sizes (38 and 70 mm in diameter) in undrained CIU (a) and drained CID (b) triaxial tests.

Figure 14 .
Figure 14.Comparison of the degradation curves G s /G 0 (a) and G t /G 0 (b) for specimens sheared at σ' 3 = const (conventional stress path) and at p' = const (vertical stress path).5.3.3.Poisson's Ratio ν and Young's Moduli E and E 50

Table 2 .
Summary of specimens for unconfined compression tests.

Table 3 .
Summary of specimens for triaxial compression tests.

Table 4 .
Average strength and stiffness parameters from unconfined compression tests.
loc-based on the readings of local on-specimen sensors; ex-based on the readings of the external sensor; * results of only 4 specimens are considered.

Table 5 .
Shear strength parameters from triaxial compression tests.

Table 6 .
Peak shear strength parameters of natural loess compacted at w ≈ w opt from other studies.